Pressure-Based System and Method for Determining Cardiac Stroke Volume 



BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 
5 [0001] This invention relates to a method for determining the stroke volume (SV) 
and, hence, any other SV-related value such as cardiac output of a human or animal 
subject, as well as to a system that implements the method. 

BACKGROUND ART 

10 [0002] Cardiac output (CO) is an important indicator not only for diagnosis of 

disease, but also for "real-time" monitoring of the condition of both human and animal 
subjects, including patients. Few hospitals are therefore without some form of 
conventional equipment to monitor cardiac output. 

[0003] One basis for almost all common CO-measurement systems is the well- 
15 known formula CO = HRSV, where SV is the stroke volume and HR is the heart rate. 
SV is usually measured in liters and HR is usually measured in beats per minute. This 
formula simply expresses that the amount of blood the heart pumps out in a minute is 
equal to the amount it pumps out on every beat (stroke) times the number of beats 
(strokes) per minute. 

20 [0004] Since HR is easy to measure using any of a wide variety of instruments, the 
calculation of CO usually depends on some technique for estimating SV. Many suitable 
techniques - both invasive and non-invasive, as well as those that combine both - are 
in use and even more have been proposed in the literature. 
[0005] One invasive way to determine cardiac output (or, equivalent^, SV) is to 

25 mount some flow-measuring device on a catheter, and then to place the catheter into 
the subject and to maneuver it so that the device is in or near the subject's heart. Some 
of such devices inject either a bolus of material or energy (usually heat) at an upstream 
position, such as in the right atrium, and determine flow based on the characteristics of 
the injected material or energy at a downstream position, such as in the pulmonary 

30 artery. Patents that disclose implementations of such invasive techniques (in particular, 
thermodilution) include: 

U.S. Patent No. 4,236,527 (Newbower et al., 2 December 1980); 



U.S. Patent No. 4,507,974 (Yelderman, 2 April 1985); 

U.S. Patent No. 5,146,414 (McKown, et al., 8 September 1992); and 

US. Patent No. 5,687,733 (McKown, et al., 18 November 1997). 



5 [0006] Still other invasive devices are based on the known Fick technique, according 
to which CO is calculated as a function of oxygenation of arterial and mixed venous 
blood. In most cases, oxygenation is sensed using right-heart catheterization. There 
have, however, also been proposals for systems that measure arterial and venous 
oxygenation non-invasively, in particular, using multiple wavelengths of light, but to date 
10 they have not been accurate enough to allow for satisfactory CO measurement on 
actual patients. 

[0007] Invasive techniques have some disadvantages, the main one of which is of 
course that catheterization of the heart is more dramatic to the patient, especially 
considering that the subjects (especially intensive care patients) on which it is 

15 performed are often already in the hospital because of some actually or potentially 
serious condition. Invasive methods also have less obvious disadvantages: Some 
techniques such as thermodilution rely on assumptions, such as uniform dispersion of 
the injected heat, that affect the accuracy of the measurements depending on how well 
they are fulfilled. Moreover, the very introduction of an instrument into the blood flow 

20 may affect the value (for example, flow rate) that the instrument measures. 

[0008] There has therefore been a long-standing need for some way of determining 
CO that is both non-invasive - or at least as minimally invasive as possible - and 
accurate. One blood characteristic that has proven particularly promising for accurately 
determining CO non-invasively is blood pressure. 

25 [0009] Most known blood-pressure-based systems rely on the so-called pulse 

contour method (PCM), which calculates as estimate of CO from characteristics of the 
beat-to-beat pressure waveform. In the PCM, "Windkessel" (German for "air chamber") 
parameters (characteristic impedance of the aorta, compliance, and total peripheral 
resistance) are used to construct a linear or non-linear, hemodynamic model of the 

30 aorta. In essence, blood flow is analogized to a flow of electrical current in a circuit in 
which an impedance is in series with a parallel-connected resistance and capacitance 
(compliance). The three required parameters of the model are usually determined 
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either empirically, through a complex calibration process, or from compiled 
"anthropometric" data, that is, data about the age, sex, height, weight, etc., of other 
patients or test subjects. U.S. Patent No. 5,400,793 (Wesseling, 28 March 1995) and 
U.S. Patent No. 5,535,753 (Petrucelli, et al., 16 July 1996) are representative of 
5 systems that rely on a Windkessel circuit model to determine CO. 

[0010] PCM-based systems can monitor CO more or less continuously, with no need 
for a catheter to be left in the patient. Indeed, some PCM systems operate using blood 
pressure measurements taken using a finger cuff. One drawback of PCM, however, is 
that it is no more accurate than the rather simple, three-parameter model from which it 

10 is derived; in general, a model of a much higher order would be needed to faithfully 
account for other phenomena, such as the complex pattern of pressure wave 
reflections due to multiple impedance mis-matches caused by, for example, arterial 
branching. Because the accuracy of the basic model is usually not good enough, many 
improvements have been proposed, with varying degrees of complexity. 

15 [001 1] The "Method and apparatus for measuring cardiac output" disclosed by 

Salvatore Romano in U.S. Published Patent Application 20020022785 A1 represents a 
different attempt to improve upon PCM techniques by estimating SV, either invasively 
or non-invasively, as a function of the ratio between the area under the entire pressure 
curve and a linear combination of various components of impedance. In attempting to 

20 account for pressure reflections, the Romano system relies not only on accurate 

estimates of inherently noisy derivatives of the pressure function, but also on a series of 
empirically determined, numerical adjustments to a mean pressure value. 
[0012] What is needed is a system and method of operation for estimating CO that is 
robust, simple, and accurate and that does not require anthropometric values or 

25 repeated calibrations. This invention meets this need. 
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SUMMARY OF THE INVENTION 
[0013] A parameter proportional to the cardiac stroke volume (SV) of a patient is 
determined by sensing an input signal that either directly indicates or is proportional to 
arterial blood pressure. The sensor used to sense the input signal may be either 
5 invasive or non-invasive. The standard deviation of the input signal is then calculated 
over a measurement interval and an estimate of SV is then calculated as a function of 
the standard deviation of the input signal. SV may be also computed as the product of 
the standard deviation and a calibration factor. In an exemplifying processing system 
that implements the method, one or more computer-executable software modules are 

10 included for carrying out the various calculations. - ~ 

[0014] Any cardiac value derived from SV may also use the invention to determine 
an SV estimate to be used for calculating the value. For example, the method 
according to the invention may be used to calculate an estimate of cardiac output (CO). 
In such an application of the invention, any known mechanism (for example, a 

1 5 hardware monitor and/or software algorithm) is used to measure the patient's heart rate 
(HR). The current cardiac output of the patient is then estimated, for example, by 
calculating the product of HR and the standard deviation and scaling the product by a 
calibration constant. 

[0015] In CO applications of the invention, the calibration constant may be 
20 determined using different techniques, both invasive and non-invasive. To calculate the 
calibration constant, a calibration cardiac output value is measured and the calibration 
constant is approved as the quotient between a calibration cardiac output estimate and 
the product of the heart rate and the standard deviation. 

[0016] The measurement interval may extend over more than one cardiac cycle, for 
25 example, to cover a time window that is multiple cardiac cycles wide. A single standard 
deviation value of the input signal may be calculated over the whole interval, or 
component standard deviation values may be calculated and then averaged (using the 
mean, median, etc.) for each of a plurality of sub-intervals to form a final composite 
standard deviation value that can be used in calculating the estimate of the cardiac 
30 stroke volume. 
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[0017] Various optimizations may be included in different embodiments of the 
invention. For example, for each of a plurality of cardiac cycles, a mean pressure value 
can be calculated and the measurement interval can then be adjusted as a function of 
change in the mean pressure value. 
5 [0018] If needed, for example, to remove the effect of potential drift in mean 

pressure over the measurement interval(s), the input signal may be high-pass filtered 
before standard deviation is calculated. 

[0019] Standard deviation may be calculated in different ways, having different 
degrees of statistical accuracy. For example, the input signal may be discretized over 
10 the measurement interval, and then a standard algorithm may be applied to determine 
an estimate of standard deviation from the sample values. As an alternative, standard 
deviation may be approximated as a function of the difference between the maximum 
and minimum values. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

[0020] Figure 1 is an illustrative example of a complex blood pressure curve over 
one beat-to-beat heart cycle 

[0021] Figure 2 is a block diagram showing the main components of a system 
according to the invention. 

20 

DETAILED DESCRIPTION 

Introduction 

[0022] In broadest terms, the invention involves the determination of a cardiac value 
such as stroke volume (SV), and/or a value derived from SV such as cardiac output 

25 (CO), using information extracted from an arterial pressure waveform, or from a 
waveform that is proportional to arterial pressure, which may be measured with an 
invasive, non-invasive, or "minimally invasive" instrument or combination of instruments. 
The organization of this description is: First, the theoretical basis of the invention is 
discussed. This is followed by an explanation of the main steps of a method to use the 

30 theory, then a description of a system that implements the method. 
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[0023] The invention may be used to advantage with any type of subject, whether 
human or animal. Because it is anticipated that the most common use of the invention 
will be on humans in a diagnostic setting, the invention is described below primarily in 
use with a "patient." This is by way of example only, however - it is intended that the 
5 term "patient" should encompass all subjects, both human and animal, regardless of 
setting. 

Theoretical Basis of the Invention 

[0024] As mentioned above, the principle formula for calculating cardiac output (CO) 
10 is CO = SV-HR, where SV is stroke volume and HR is heart rate. Given HR, the 
problem then remains how to determine SV. Based on the observation that the 
pulsatility of a pressure waveform is created by the cardiac stroke volume into the 
arterial tree, inventors have discovered that one particularly elegant solution to the 
problem is to estimate SV to be proportional to the standard deviation of the arterial 
15 pressure waveform P(t), or of some other signal that itself is proportional to P(t). Thus, 
according to this aspect of the invention, 

SV = K-Gp 
from which follows that 

CO = KopHR 

20 where K is a constant and a P is the standard deviation of P(t) (or of some other signal 
proportional to P(t)) taken over some time interval. Other quickly computed functions of 
P(t) that provide values proportional to SV are also discussed below. 

Steps for Determining CO 

25 [0025] As mentioned above, the invention may be used to determine SV, and hence 
any other cardiac value that is a function of SV. It is anticipated that the most common 
use of the invention will be to determine CO, so the exemplary embodiments of the 
invention are described below primarily with respect to this application. Since the 
invention generates an estimate of SV, however, those skilled in the art of medical 

30 systems will know how to use the estimated SV value to derive still other cardiac 
values. 
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[0026] The first step in the method for measuring CO is to acquire a representation 
of either the arterial pressure waveform P(t) or of some waveform that is proportional to 
P(t). This may be a direct measurement of arterial pressure, or a measurement of 
some other parameter that is proportional to arterial pressure. For example, a 
5 plethysmographic technique, for example, using a finger cuff, will produce a signal that 
is proportional to arterial blood pressure; this signal can then be scaled to convert it to 
an estimate of blood pressure. For the sake of succinctness, the method for 
determining SV is discussed below with reference to the arterial pressure waveform 
P(t), but one should keep in mind that the disclosed method steps may easily be 

10 adjusted to use another, proportional waveform, simply by applying the proper scaling, 
such as the scaling provided by a calibration factor K described below. 
[0027] Figure 1 illustrates an example of the waveform P(t) of arterial pressure taken 
over a single heart cycle, here, from the point of diastolic pressure P dia at time t^, 
through the time t sys of systolic pressure P sys , to a time t^, at which the blood pressure 

15 once again reaches P dia . 

[0028] According to the invention, P(t), or any signal that is proportional to P(t), may 
be measured at any point in the arterial tree, either invasively or non-invasively. If 
invasive instruments are used, in particular, catheter-mounted pressure transducers, 
then any artery may be used as a measurement point. Placement of non-invasive 

20 transducers will typically be dictated by the instruments themselves - the placement of 
finger cuffs, upper arm pressure cuffs, and earlobe clamps should be obvious. 
Regardless of the instrument, it will ultimately produce, or cause to be produced, an 
electric signal corresponding (for example, proportional) to P(t). 
[0029] As is well known, analog signals such as P(t) can be digitized into a 

25 sequence of digital values using any standard analog-to-digital converter (ADC). In 
other words, P(t), tO < t < tf, can be converted, using known methods and circuitry, into 
the digital form P(k), k = 0, (n-1 ), where tO and tf are initial and final times, respectively, 
of the measurement interval and n is the number of samples of P(t) to be included in 
the calculations, distributed usually evenly over the measurement interval. 

30 [0030] The next step of the procedure is to calculate the standard deviation of P(k), 
that is, c P , or some value that has a known relationship to the standard deviation (see 
below). One way to calculate a P is to use the well-known algorithm for calculating the 
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standard deviation of a discretized function. Thus: 

4 -- ~ £ [P(k) -P a ] 2 (Equation 1 ) 

where P avg is the mean pressure value, that is: 

Pavg=-ZP(k) (Equation 2) 

n k=i 

Of course, to get a P the system simply takes the square root of a P 2 . 

[0031] The analog measurement interval, that is, the time window [to, tf], and thus 
the discrete sampling interval k=1,n, over which a CO estimate is calculated should be 
small enough so that it does not encompass substantial shifts in the mean pressure 
P avg . Also, one could filter out low frequency variations such as respiration using a high 
pass filter, which would also help remove the effect of any drift in mean arterial pressure 
during the time window. For the sake of providing more stable and reliable readings, 
however, is it best to let the time window for each CO estimate extend longer than one 
cardiac cycle. Preferably, the measurement interval (time window) should be a plurality 
of cardiac cycles, that is, beginning and ending at the same point in different cardiac 
cycles; this ensures that the mean pressure value used in the calculations of a will use 
a mean pressure value P avg that is not skewed because of incomplete measurement of 
a cycle. Since the information is available to any embodiment of the invention that uses 
Pavg (for example, in the calculation of a P ), it may be useful in some implementations to 
always display the current mean blood pressure and/or trend to the user. 
[0032] Larger sampling windows have the advantage that the effect of perturbations 
such as those caused by reflections will usually be reduced, since they will be tend to 
"cancel out" in the calculations of mean pressure and standard deviation. An 
appropriate time window can be determined using normal experimental and clinical 
methods. Note that it would be possible for the time window to coincide with a single 
heart cycle, in which case mean pressure shift will not be of concern. 
[0033] As a check, the system according to the invention could also, as a separate 
background operation, compute the mean pressure over each cardiac cycle. If the 
mean cycle-to-cycle pressure shows any absolute or proportional drift greater than 
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some threshold value, a warning signal could be generated such that the currently 
computed CO estimate may be considered less reliable or discarded altogether. 
[0034] It would be also possible to adjust the time window [tO, tf] according to drift in 
P avg . For example, if P avg over a given time window differs absolutely or proportionately 
5 by more than a threshold amount from the P avg of the previous time window, then the 
time window could be reduced; stability of P avg could then be used to indicate that the 
time window can be expanded. The time window could also be expanded and 
contracted based on noise sources, or on a measure of SNR or of SNR variation. In a 
preferred embodiment, limits are placed on how much the time window is allowed to 

10 expand or contract; and an indicatiorf of the time interval is displayed to the user. 
[0035] It is not necessary for the time window to start at any particular point in the 
cardiac cycle. Thus, to need not be the same as t^, although this may be a convenient 
choice in many implementations. This means that the beginning and end of each 
measurement interval (that is, tO and tf), each yielding a CO estimate, may be triggered 

15 on almost any characteristic of the cardiac cycle, such as at times t^ or t sys , or on non- 
pressure characteristics such as R waves, etc. 

[0036] Rather than calculate a single a P value from a multi-cycle measurement as in 
Equation 1 it would also be possible to calculate several ct p , for example, one for each 
of a plurality of intervals, and then to average them (by taking the mean, median, etc.) 
20 to compute a composite a P for use in the formulas. The inventors have, moreover, 
discovered other alternatives for computing a pulsatility variable similar to the standard 
deviation a P of the arterial blood pressure. 

[0037] The inventors have observed that other values may be derived from the 
pressure waveform P(t) that either provide an approximation of o P , or that also are 

25 proportional to SV, or both. As one example, the inventors have observed that the 
difference between the maximum and minimum measured pressures, taken over the 
time window, is a pulsatility measurement that may be substituted for direct calculation 
of a P using the standard formulas given above as Equations 1 and 2. Let max[P(k)] and 
min[P(k)] be the maximum and minimum values, respectively, of the sampled pressure 

30 over the measurement interval. The standard deviation is approximately equal to one- 
third times the difference of these two values: 
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a P * {max[P(k)] - min[P(k)]} /3 



[0038] Although probably less accurate than calculation of a P using Equations 1 and 
2, this "rough" a P approximation has the advantage of simplicity, requiring no sampling 
5 of P(t) at all. Indeed, given an input signal indicating heart rate (HR), a system to 

compute (max[P(k)] - min[P(k)]} /3 and, from it, SV and/or CO (or some other function of 
SV) could be implemented completely in hardware, even all-analog circuitry, using 
known circuit design techniques. This would allow development of very inexpensive, 
easily manufactured and physically robust CO monitors for use in areas or applications 
10 that have only minimal facilities and resources. 

[0039] Of course, it is not necessary to have a separate calculation relating 
max[P(k)] and min[P(k)] to a P , and then to use a P to calculate SV. This is described 
here by way of illustration only. Rather, given max[P(k)] and min[P(k)], SV can be 
estimated directly as: 

15 

SV = k {max[P(k)] - min[P(k)]} 

where k = K/3. (Of course, K can simply be adjusted to account for the factor 1/3.) 

20 [0040] As another alternative, the inventors have observed that the maximum or 
absolute value of the minimum of the first derivative of the P(t) with respect to time is 
generally proportional to a P and to SV. Thus: 

■ rdP(t)^ 
mm — — 
I dt J 

25 [0041] It would also be possible to use the average or these first derivatives instead 
of using only the one or the other. Given P(k), the derivatives may be determined using 
any known numerical method; note that the points of interest on the pressure waveform 
are the points of inflection, that is, the points at which the second time derivative of P(t) 
is zero. The time interval over which these derivatives is evaluated may be the entire 

30 cardiac cycle. It will generally suffice, however, to evaluate P(t) between the beginning 
of the cardiac and the first dicrotic point, shown as P djcrotic in Figure 1 , since the 

10 



SV=Kmax 



f dP(t) 
I dt 



or SV=K 



maximum positive slope will usually occur about halfway between the diastolic and 
systolic points, that is, P dia and P sys and the greatest negative slope will generally occur 
about halfway between the systolic and first dicrotic points, that is, P sys and P dicrotic . 
Examining only these portions of P(t) will eliminate the possibility that spurious values 
5 will be used from after the time of P dicrotic . 

[0042] As yet another alternative, a standard software or hardware module could be 
used to compute the Fourier transform of the measured pressure signal P(t) over each 
cardiac cycle, or over a multiple of cycles. The quantity defined by the magnitude of the 
Fourier transform component, H1, at the primary frequency component, i.e., the 
10 frequency corresponding to the "heart rate," divided by the mean arterial pressure P avg , 
that is, H1/P avg , will be proportional to SV. Instead of H1, the magnitude of the Fourier 
component H2 corresponding to twice the heart rate, that is, the first harmonic 
frequency component, cold be used instead; thus, H2/P avg will also be proportional to 
SV. 

1 5 [0043] In order to calculate CO, the heart rate HR (or some signal from which HR 
can be derived) is needed. Any of the many known instruments for measuring HR may 
be used. If the beginning and end times for each P(t) interval are triggered by an 
electrocardiogram signal, for example, then the same signal may be used to calculate 
HR. The measured pressure wave P(t) (in practice, P(k)) may itself be used to derive 

20 HR, for example, using standard Fast Fourier transformation or derivative analysis. 
[0044] Before finally arriving at a value for CO, it is also necessary to determine a 
value for the calibration constant K. One way to do this is as any pre-determined 
function of P(t); thus, K = K(P(t)). In this case no independent CO technique is 
necessary. 

25 

[0045] Another way to do this is to use any known, independent CO technique to 
determine this relationship, whether invasive, for example, thermodilution, or non- 
invasive, for example, trans-esophageal echocardiography (TEE) or bio-impedance 
measurement. The invention provides continuous trending of CO between intermittent 
30 measurements such as TO or TEE. Using the chosen independent method, a value 
CO^, is determined, so that K will be: 
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K = CO ca /(V.HR) 



where V is the chosen value proportional to SV, for example: 
V=a P; or 

V= max[P(k)] - min[P(k)] ; or 

V= maximum or absolute value of the minimum of the first derivative of the P(t); or 
V=H1/P avg orH2/P avg 



[0046] Even if an invasive technique such as catheterization is used to determine K, 

10 it will usually not be necessary to leave the catheter in the patient during the 

subsequent CO-monitoring session. Moreover, even when using catheter-based 
calibration technique to determine K, it is not necessary according to the invention for 
the measurement to be taken in or near the heart; rather, the calibration measurement 
could be made in the femoral artery. As such, even where an invasive technique is 

15 used to determine the calibration constant K, the invention as a whole is still minimally 
invasive in that any catheterization may be peripheral and temporary. 
[0047] As is mentioned above, rather than measure arterial blood pressure directly, 
any other input signal may be used that is proportional to blood pressure. This means 
that calibration may be done at any or all of several points in the calculations. For 

20 example, if some signal other than arterial blood pressure itself is used as input, then it 
may be calibrated to blood pressure before its values are used to calculate standard 
deviation, or afterwards, in which case either the resulting standard deviation value can 
be scaled, or the resulting SV value can be calibrated (for example, by setting K 
properly), or some final function of SV (such as CO) can be scaled. In short, the fact 

25 that the invention may in some cases use a different input signal than a direct 

measurement of arterial blood pressure does not limit its ability to generate an accurate 
SV estimate. 



System Components 

30 [0048] Figure 2 shows the main components of a system that implements the 
method described above for sensing pressure and calculating CO. As is mentioned 
above, pressure, or some other input signal proportional to pressure, may be sensed in 
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either or, indeed, both, of two ways: invasively and non-invasively. Simply because it is 
anticipated to be the most common implementation of the invention, the system is 
described as measuring arterial blood pressure (as opposed to some other input signal 
that is converted to pressure) and generating a CO estimate (as opposed to just SV). 
5 Changes to the illustrated system to accommodate other design choices will be obvious 
to those skilled in the art of medical devices. 

[0049] Figure 2 shows both types of pressure sensing for the sake of conciseness; 
in most practical applications of the invention, either one or several variations will 
typically be implemented. In invasive applications of the invention, a conventional 

10 pressure sensor 100 is mounted on a catheter 110, which is inserted in an artery 120 of 
a portion 1 30 of the body of a human or animal patient. Such artery could be an 
ascending aorta, or pulmonary artery, or, in order to reduce the level of invasiveness, 
the artery 120 could be peripheral, such as the femoral, radial or brachial artery. In the 
non-invasive applications of the invention, a conventional pressure sensor 200, such as 

15 a photo-plethysmographic blood pressure probe, is mounted externally in any 

conventional manner, for example using a cuff around a finger 230 or a transducer 
mounted on the wrist of the patient. Figure 2 schematically shows both types. 
[0050] The signals from the sensors 100, 200 are passed via any known connectors 
as inputs to a processing system 300, which includes one or more processors and other 

20 supporting hardware and system software (not shown) usually included to process 
signals and execute code. The invention may be implemented using a modified, 
standard, personal computer, or it may be incorporated into a larger, specialized 
monitoring system. In this invention, the processing system 300 also may include, or is 
connected to, conditioning circuitry 302 which performs such normal signal processing 

25 tasks as amplification, filtering, ranging, etc., as needed, as well as the optional high 
pass filtering mentioned above. The conditioned, sensed input pressure signal P(t) is 
then converted to digital form by a conventional analog-to-digital converter ADC 304. 
As is well understood, the sampling frequency of the ADC 304 should be chosen with 
regard to the Nyquist criterion so as to avoid aliasing of the pressure signal; this 

30 procedure is very well known in the art of digital signal processing. The output from the 
ADC 304 will be the discrete pressure signal P(k), whose values may be stored in 
conventional memory circuitry (not shown). 
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[0051] The values P(k) are passed to (usually, accessed from memory by) to an SV- 
calculation module 306, which is a software component comprising processor- 
executable code for calculating whichever value V is used to determine SV as 
explained above. For example, where a P is calculated directly, the SV-calculation 
5 module 306 will evaluate Equations 1 and 2 above, or equivalent expressions. The 
calculation module 306 preferably also selects the time window [to, tf] over which each 
CO estimate is generated. This may be done as simply as choosing which and how 
many of the stored, consecutive, discretized P(t) values P(k) are used in each 
calculation, which is the same as selecting n in the range k = 1, n,_ . 

10"" [0052] The computer value of SV is then passed to a subsequent CO-calculation 
module 308, which similarly comprises executable code for evaluate the expression 
CO = K-Gp-HR. Of course, modules 306 and 308 may be combined into a single 
software component; they are shown separately for the sake of clarity. Of course, the 
CO-calculation module 308 requires values for HR and K as well: 

15 [0053] The patient's current heart rate HR is either calculated from the measured 
pressure curve P(k) by a corresponding software module 310 (for example, using 
Fourier or derivative analysis) or is otherwise measured with any conventional hardware 
device. 

[0054] The value K will normally be input to module 308 automatically, but can be 
20 entered by the operator via a conventional input device 400 such as a keyboard, 
mouse, etc. This will usually be the same input device(s) used by the processing 
system 300 for other purposes such as entering data identifying the patient and the 
specifications of the monitoring session. 

[0055] The module 308 calculates and estimates CO for each chosen measurement 
25 interval. Each estimated CO value is preferably output to any conventional display or 
printing device 500 for the user to view and monitor. As with the input device 400, the 
display 500 will typically be the same as is used by the processing system for other 
purposes. 

[0056] The invention further relates to a computer program loadable in a computer 
30 unit or the processing system 300 in order to execute the method of the invention. 
Moreover, the various software modules 306, 308, 310, and 312 used to perform the 
various calculations and perform related method steps according to the invention may 



14 



also be stored as computer-executable instructions on a computer-readable medium in 
order to allow the invention to be loaded into and executed by different processing 
systems. 

Other outputs 

[0057] The invention is described above in the context of calculating estimates of 
CO. This is the use of invention that the inventor assumes will be most common, but 
the invention is not inherently limited to such use. In essence, the invention provides a 
novel way to calculate stroke volume SV, or any parameter that is a function of (for 
example, proportional to) SV, not just CO. Consequently, the advantages of the 
invention will apply to the calculation of any value derived from SV. For example, the 
end diastolic volume (EDV) and the ejection fraction (EF) are related as EF=SV/EDV, 
which expresses the intuitive relationship that the pumping efficiency (EF) of the heart is 
the ratio between how much blood the heart pumps out on every beat (contraction) and 
how much blood is in the heart chamber just before the beat. Inversely, EDV=SV/EF. 
If an estimate for either EDV or EF is determined in some other known manner, then 
the invention could be used to provide SV, and thus an estimate of the other of EDV or 
EF. 
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